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Abstract 

The interaction of screw dislocations with an applied stress is studied using atomistic 
simulations in conjunction with a continuum treatment of the role played by the far field 
boundary condition. A finite cell of atoms is used to consider the response of dislocations 
to an applied stress and this introduces an additional force on the dislocation due to the 
presence of the boundary. Continuum mechanics is used to calculate the boundary force 
which is subsequently accounted for in the equilibrium condition for the dislocation. Using 
this formulation, the lattice resistance curve and the associated Peierls stress are calculated 
for screw dislocations in several close packed metals. As a concrete example of the boundary 
force method, we compute the bow out of a pinned screw dislocation; the line-tension of 
the dislocation is calculated from the results of the atomistic simulations using a variational 
principle that explicitly accounts for the boundary force. 

1 Introduction 

Atomistic simulations with carefully chosen potentials have been shown to account for the struc- 
ture and energetics of a range of extended defects (Vitek 1992). Despite these successes, one of 
the main drawbacks of atomistic simulations is that they are limited to small simulation cells, 
necessitating an approximate treatment of boundary conditions. A popular method used in the 
simulation of dislocations is to use a finite sized simulation cell and to adopt the continuum 
linear elastic solution for the defect of interest as the initial trial solution for the atomistic energy 
minimisation with the additional constraint that the atoms on the boundary of the cell be fixed at 
the positions dictated by the elastic solution. The rationale behind this approach is that at large 
distances from the core of the defect, continuum linear elasticity provides an accurate description 
of the displacement, strain and stress fields. This method has been applied with much success to 
the determination of the static core structure of dislocations (Vitek 1992). Upon application of 
an applied stress, the defect of interest will adopt a new equilibrium configuration. In particular, 
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the dislocation core moves relative to the axis of symmetry of the simulation cell. This motion 
is influenced by the boundary, i.e., the defect encounters an additional force that depends on its 
position relative to the boundary of the cell. Failure to account for this boundary effect in the 
analysis of the interaction of the defect with an applied stress leads to results that are simulation 
cell-size dependent making the quantitative determination of mechanical parameters such as the 
line tension and the Peierls stress of the dislocation difficult. Indeed, Basinski, Duesbery and 
Taylor (1971) observed such effects in their study of the interaction of screw dislocations with 
an applied stress in a sodium lattice. Although the importance of the boundary stress has been 
realised, a systematic analysis of its effects has been wanting. 

In this paper we develop a method that explicitly accounts for the boundary force using linear 
elasticity, for the case of screw dislocations. It is shown that the boundary force scales as the 
inverse square of the size of the simulation cell. This formulation is used in conjunction with 
atomistic simulations to obtain the lattice resistance function for screw dislocations in several 
close packed metals. In addition, the bow out of a pinned screw dislocation under the influence 
of an applied stress is studied using a variational principle that incorporates the boundary force, 
demonstrating that in this case the synthesis of atomistic and continuum methods yields a 
consistent description of bow out. One outcome of this calculation is that the line tension of the 
dislocation can be computed directly from atomistics. In the following section we present the 
derivation of the boundary force after discussing its origin and significance. Section 3 contains 
the results of the lattice resistance calculations, while the atomistic study of dislocation bow out 
may be found in Section 4. We conclude with a discussion of the significance of our findings. 

2 Simulation methodology and the boundary force 

Our objective is the simulation of the interaction of an infinite straight screw dislocation, in an 
infinite crystal, with an applied stress. To model this problem using a finite simulation cell we 
choose a cylindrical cell of atoms with axis along the ^-direction (cf. fig. 1). The X\Xi plane 
is the slip plane, the ^-direction is the line direction as well as the Burgers vector direction 
of the dislocation. A screw dislocation is placed at the center of this cell at O by displacing 
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the atoms from their perfect lattice sites according to the continuum linear elastic solution 
for a Volterra dislocation and this configuration is used as the trial solution for the atomistic 
energy minimisation. The minimisation is performed keeping the atoms in region F fixed at the 
positions given by the linear elastic solution; only the atoms in region M [dynamic region, a 
cylinder of radius R) enclosed by F are allowed to move. The relaxed core of the dislocation 
is thus obtained. The interaction of the dislocation with a homogeneous applied shear stress is 
simulated by displacing the atoms in both F and M by the displacement fields corresponding to 
the homogeneous strain that is equivalent to the stress of interest. Again, the atomistic energy 
minimisation step is performed keeping the atoms in the region F fixed and allowing the atoms 
in the region M to move. On performing the above procedure, the dislocation moves through 
the crystal to a new position P, which is determined by the equilibrium condition 

Tapp + T b {d) +T L {d) = 0, (1) 

where d = OP as shown in fig. 1, r a , pv is the resolved shear stress which effects a Peach-Koehler 
force T app b in the Xi-direction (b is the magnitude of the Burgers vector), tl is the lattice resistance 
function (Kocks, Argon and Ashby 1975) which is periodic in d and u is the boundary stress 
which is a monotonically increasing (in magnitude) function of d. 

We now explore the origin and significance of the boundary stress term. Upon application 
of a homogeneous strain, the net displacement of the atoms before relaxation in region F and 
on the external boundary dM e of M is equal to the sum of the elastic displacements due to 
the dislocation at O (uo) and the displacements due to the homogeneous strain as dictated by 
linear elasticity. After relaxation, the dislocation moves to point P; if the dislocation was set 
in an infinite crystal, the displacements of the atoms in region F would have been equal to 
the sum of the elastic displacements due to the presence of the dislocation at P (lip) and the 
displacements due to the homogeneous strain. Since the atoms in the region F are held fixed 
during the relaxation, the boundary condition after relaxation is not consistent with the fields of 
the dislocation at P. This inconsistency causes a spurious increase in the total energy stored in 
the system, i.e., there is an additional energy cost in moving the dislocation from O to P. This 
additional energy produces an energetic stress on the dislocation (in the sense of Eshelby 1951) 
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which we call the boundary stress. This boundary stress depends on the radius R of the dynamic 
region, making the quantities calculated from the results of the simulations, cell-size dependent. 
The boundary stress is particularly important in the simulations of dislocations in crystals with 
low lattice resistance such as those with the fee and hep structures. In materials with high lattice 
resistance such as bec metals this term may be neglected as the lattice resistance term will be 
several orders of magnitude larger than the boundary term. 

In order to derive the boundary stress, we assume that the total energy S is built up of 
an elastic energy E and a misfit energy M, i.e., S = E + M. The misfit energy accounts 
for the energy required to maintain the slip distribution due to the presence of the dislocation, 
and the elastic energy is due to the elastic straining of the crystal and may be thought of as 
arising from a continuous distribution of infinitesimal dislocations. Thus the net change in the 
energy (AS) on moving the dislocation to the point P is then given by AE(d) + AM.(d). In 
the absence of the inconsistent boundary condition, the change in elastic energy would vanish 
( in accordance with a Peierls-Nabarro type model which states that the elastic energy remains 
unchanged as the dislocation moves through the crystal). We assume that the presence of the 
inconsistent boundary does not effect the misfit energy, and thus the additional energy due to 
the inconsistent boundary is entirely elastic. In addition, we assume that the change in elastic 
energy AE computed using a Volterra model of the dislocation gives a good approximation to 
this energy change. These assumptions are valid when the size of the simulation cell (R) is large 
compared to the core width of the dislocation. Under these assumptions, the change in elastic 
energy AE may be computed and the boundary stress is given by 

1 dAE 

^ = -s-sp (2 > 

and the lattice resistance function t l may be obtained as 

TL{d) = -l-dT- (3) 

In the following discussion, the external boundary of the dynamic region is denoted by dM e , the 



slip surface AO is denoted by dM s and dM s > denotes the slip surface AP that corresponds to the 
dislocation at P. All the linear elastic fields associated with the presence of the dislocation at 
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O are denoted by a subscript O (e.g. u ) while a subscript P is used to denote fields associated 
with the dislocation at P. In both cases, the fields correspond to an unconstrained dislocation in 
an infinite crystal. The dislocation is displaced from O to P, while keeping the displacements on 
dM e fixed at u G (a constrained motion of the dislocation). This corresponds to the process that 
occurs in the atomistic simulation when the dislocation moves under the influence of the applied 
stress. It is clear that AE is the difference in the elastic energy E 2 of the system when the 
dislocation is at P after the constrained motion and the elastic energy E\ when the dislocation 
is at O. The energy E\ is evaluated as 

E x = \\ t - hojdS +\j t - u Q dS, (4) 

2 JdM 3 2 JdMe 

where t Q is the traction vector on the relevant boundary due to the stress field of the dislocation 
at O ( to = (To • n where cr is the stress tensor and n is the normal to the boundary/surface), 
and [ ] represents the jump in the displacement across the slip surface (i.e., the Burgers vector). 
The total energy E 2 of the configuration with the dislocation at P (after the constrained motion) 
may be evaluated as follows. This configuration is taken as the superposition of the fields u P 
due to the presence of the dislocation at P as though in an infinite crystal, and Au which is an 
equilibrium displacement field with 

Au = u - up, Vx G dM e (5) 

and with the condition that [Au] = on dM s i. It must be noted that Au at a point in the 
interior of M is not equal to uo — up. The energy in the displaced configuration is 

E 2 = -[ tp • [iipldS + - / tp • updS + [ t P ■ AudS + - / At • AudS, (6) 

2 JdM s , 2 JdM e JdM e 2 JdMe 

where tp is the traction associated with the stress field of the dislocation at P (as though in 
an infinite body), At is the traction associated with the fields Au. The additional energy AE 
stored in the system due to the boundary condition on dM e may be expressed as 

AE = E 2 -Ei = AE a + AE b + AE C , (7) 

with the three terms being given by 

AE a = \f t P -{u P \dS-\ ( to -holdS, (8) 

2 JdM, 2 JdM s 



AE b = f t P • (uo - ^u P )dS - \ I t - u Q dS, (9) 

JdM e 2 2 J9M e 



AE C 



= - I At • AudS. (10) 

2 .W V ' 



We choose to partition the energy change in this fashion for convenience as it will be shown 
below how these three contributions may each be evaluated separately 

In preparation for evaluating the energy difference, we record the isotropic, linear elastic 
displacement and stress fields for a Volterra screw dislocation (Hirth and Lothe 1968) at a point 
Q with coordinates (xq,0) 

(u 2 )q(x u x 3 ) = u Q (x 1 ,x 3 ) = -^-arctan ( — — — J. (11) 

Z7C \Xl — Xq J 

The point Q is a generic point on the slip plane; Q may, for example, be taken to be O or P. 
The components of the stress tensor <tq are given by 

where = arctan(x 3 /(a; 1 — xq)), r' = ^ (x~i — xq) 2 + x 3 and /x is the shear modulus. All 
other displacement and stress components not symmetric to those given above vanish. In the 
developments that follow, we freely switch between cartesian (x±, x 3 ) and polar coordinates (r, tp). 
The notation is simplified by adopting the convention that 

oq = t Q ■ e 2 = (ctq • n) • e 2 = (a^QUi + (a 32 )Qn 3 , (13) 

th 

where n is the normal to the pertinent boundary and denotes the basis vector along the i 
coordinate axis. 

First we compute AE a which requires the evaluation of surface integrals over the domains 
dM s and dM s i for which n = e 3 . Taking r Q to be the elastic core cut-off radius, we see that 

\[ t P ■ {upjdS = ^- In - + d , and J f t • [uo]dS = ^- In — . 

2 JdM s , 4:7V r 2 J dM 3 4:71 r Q 

It is now immediate that the contribution of these terms to the energy difference is 

AE a = ln(l + a), (14) 
4ir 



with a = d/R. 

To evaluate AE b we note that the unit normal to the external boundary dM e is given by 
n(-0) = cos^e! + sin-?/>e3. Using eqs. [12] and [H| it follows that to ■ e 2 = <j q = on dM e and 
hence from eq. (Qj 

AE b = [ t P ■ (u - }-up)dS = [ a P (u - \u P )dS. (15) 

JdM e 2 JdM e 2 

Since in this case xp = d, one obtains that 

fibd sin ip 



a P (R,i>) 



2tt r' 

for points with coordinates {R,ip) on 9M e noting that = arctan (sin?/'/ (cos ^ — a)) and r' 
y/R' 2 + d 2 — 2Rd cos ■0. Substituting all relevant quantities into eq. ( |I5|) , results in 

/i6 2 



(16) 



4vr 2 



al\(a) 



(17) 



where 



Ji(a) 



,1 f smtb 

%p arctan 

2 \ cos-0 — a 



sin-0 



V2 



and p' 2 = (1 + a 2 — 2a cos 0). 

We now turn to the evaluation of AE C . From eq. (|10|) we have that 

AE C = - ( At • AxidS = - [ AaAudS 

2 JdM e 2 JdMe. 



(19) 



with Au = Aue2 and (Act ■ n) ■ n = Act. The evaluation of AE C , requires both knowledge of 
An on dM e , which is provided by boundary conditions, and of Act. Since Au is an equilibrium 
field, it satisfies the elastic equilibrium equation (Navier's equation) in M, 



V 2 (Am) 



with the associated boundary condition that on dM e 
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Au = arctan — — arctan 

27T V \Xi J \Xi — d 



(20) 



(21) 



This Dirichlet problem for Am in the interior of M may be solved using the Poisson integral 
formula (John 1982), i.e., 



Au(r, xjj) 



4vr 2 



9 — arctan 



sin 6 



R 2 -r 2 



cos 9 — a J J \R 2 + r 2 — 2rR cos (0 — 



dO 



(22) 
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where (r, if>) are the polar coordinates of an interior point of M. To evaluate Aa on dM e we 
evaluate the strain field corresponding to the displacement field given by eq. (^) and use the 
elastic constitutive equations. This results in 

pb 



Aa 



2ix 2 R 



j(V0 



(23) 



with J (if}) defined by 



j(V0 



= lim 
= lim 



2vr 2 «9Au(p, V) 



dp 

',-axctan 



sin 



cos 6* — a 



' 2p - (p 2 + 1) cos (if) - 6) 
(l + p 2 -2pcos (if)-9)f 



where p is the normalised value of r, i.e., p = r/R. Substituting eqs. |21] and |23| into eq. 
find 



AE C 



ptf 



h(a) 



(24) 



we 



(25) 



with 



Ha) 



J(ip) if) — arctan 



sin^ 



dip. 



(26) 



^cosip — a / 

Having evaluated the various terms that contribute to the elastic energy as a result of the 
inconsistent boundary condition, we are now in a position to evaluate the boundary stress. 



Collecting terms AE a , AE and AE C from eqs. [14], [17] and |25] one obtains that 

ph 2 



AE 



ATI 2 



Z(a) 



(27) 



where Z(a) is given by 



Z(a) = 7iTn (1 + a) — al\(a) /2(a). 



(21 



It is remarkable that AE does not depend explicitly on d and R, but depends only on their ratio 
a. Also, the result is independent of the elastic core cutoff radius r a (on the assumption that it 
remains fixed). 

Although the above discussion was restricted to isotropic linear elastic solids, our method 
may be applied to the anisotropic case as well. We indicate the changes that are required to 
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accomplish this in the case of hep crystals with (1120) type screw dislocations. In this case the 
simulation cell is chosen to be an elliptic cylinder, i.e., the boundary dM e satisfies the equation 

x\ + |f = R 2 , (29) 



where k = yC^/C^ and Cy have their usual meanings (Hirth and Lothe 1968). The terms 
AE a and AE b can be evaluated straightforwardly by following the steps of the previous section. 
To evaluate AE C we note that Au satisfies 

with boundary conditions as in eqs. |5|. This equation may be solved after a coordinate transfor- 
mation (x\, £3) — > 1/3) where y% = x\ and y% = x^/k; Au is harmonic in the new coordinates 
and dM e transforms to a circle of radius R, and thus Am may be calculated using the Poisson 
integral formula. The traction Aa is calculated using the elastic constitutive equations and AE C 
is evaluated as before. The final expression for AE may be written in this case as 

AE = ^Z(a,k) (31) 

where in this case the dimensionless energy change Z depends on both a and the ratio k which 
is a measure of the elastic anisotropy. The case of isotropic solids is obtained from eq. (|3T| ) by 
using k — 1.0 and C44 = fi. 

The boundary stress may now be obtained using eq. (Q) as 

C u b dZ(a,k) C u b , M , , 

Tb{a) = ~^R -^T~ = " } ' (32) 

where E(a, k) is the derivative of Z with respect to a. The dimensionless stress F is plotted in 
fig. 2. For small values of a, Z admits a quadratic approximation of the form 

Z{a,k) = ^A{k)a 2 , (33) 

which demonstrates that when a C 1 (i.e., d R) the boundary stress is linear in the dimen- 
sionless parameter a. This linear nature of F in the indicated regime is evident from fig. 2. The 
constant A, which depends on the value of k, is the boundary stress coefficient and the boundary 
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stress u becomes 

C M bA(k) 

The values of the boundary stress coefficient are computed to be 6.28 for the isotropic case and 
6.19 for the case of titanium (HCP). It is clear from eq. ([34] ) that the boundary force is directly 
proportional to d for a given R and for a given d depends on the inverse square of the size of the 
simulation region M. Assuming that a typical Burgers vector is of the order of 5.0A and that 
the radius of the simulation cell is about 50. OA one obtains a boundary stress of about 0.00025/i, 
when d is taken to be one Burgers vector. Thus the boundary effects become highly significant 
in the cases of materials with low lattice resistance (fee and hep metals). On the other hand, 
in the case of BCC screw dislocations, the boundary term may be neglected in comparison with 
the Peierls stress which is two orders of magnitude larger. 

In light of our development of the boundary stress, we now proceed to derive the equilibrium 
condition (eq. (HD). We assume that the dislocation has moved to the point P under the influence 
of the applied stress r app and is in equilibrium. The principle of virtual work states that if the 
dislocation is given an arbitrary virtual displacement S from P, the net external work done will 
equal the change in the net energy stored (per unit thickness along the ^-direction), i.e., 

dAE dAM 

Tappbd = —^j-o H ^—d . (35) 



dd dd 

external work 



change in energy 
Eq. |35| can be rewritten as 

{T app + n{d) + r L {d))8 = 0. (36) 

Making use of eqs. and [3[ and noting that S is arbitrary, the equilibrium condition (eq. ([I])) is 
obtained. By virtue of the continuum arguments outlined here, we are now in a position to use 
the boundary force in the context of explicit atomistic simulation. 

It may be noted that, in the isotropic case, the boundary force can derived by placing an image 
dislocation at a point (R 2 /d,0) on the X\ axis. The resulting fields obtained by superposition 
will then satisfy the condition eq. fl5|) and the boundary force is now equal to the force exerted 
by the image dislocation on the dislocation at P. The advantage in using the energetic method 
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of obtaining the boundary force is the generalization to the anisotropic case, and possibly for 
other defects such as the edge dislocation. 

3 Lattice resistance functions and Peierls stress 

In this section we present the results of the calculation of the lattice resistance function for screw 
dislocations in several close packed metals. As mentioned above, the form of the lattice resistance 
function was discussed by Kocks, Argon and Ashby (1975). Although atomistic calculations have 
been used previously to compute the Peierls stress which corresponds to a single point on the 
lattice resistance curve, it is one of the aims of the present work to compute the entire lattice 
resistance function from atomistic simulation. This is achieved using the following method. The 
position of the dislocation (d) is monitored as a function of the applied stress. The applied stress 
is sufficiently small so that the displacement d is small compared with R; thus the linear boundary 
stress approximation may be applied (eq. (|34D). From the equilibrium condition (eq. (|l|)) the 
lattice resistance function can be obtained as 



The motion of the dislocation is monitored relative to an initial equilibrium position labeled 
as point O which is chosen to correspond to a Peierls valley. To obtain the position d of the 
dislocation the following procedure is adopted. The slip distribution on the slip plane (xix 2 ) i.e., 
distribution of [7/2] as a function of x\ is obtained using the method given by Miller and Phillips 
(1996). The position d of the dislocation is taken to be that value of x\ at which [iz 2 ] equals 
one half the Burgers vector (6/2). For this procedure to be meaningful it is essential that as the 
dislocation moves through the crystal, the slip distribution moves in an essentially self similar 
fashion, which we find is a good approximation in the cases considered here. 

To obtain the relationship between the applied stress r app and d, the stress is applied in 
small increments and the configuration is relaxed at every step as explained before. Once the 
applied stress reaches a desired maximum value (this is chosen so that the resulting d/R value 
is not overly large), the stress is decremented until it reaches zero, i.e., the direction of loading 
is reversed. Fig. 3 shows the net resistance that the dislocation would have to overcome in its 




(37) 



11 



journey through the crystal (i.e., + tl). The expected r app vs. d curve is OABCDEFG as 
shown in fig. 3 in the loading branch and GHIJKLMO during the unloading process. Thus if 
the applied stress is sufficiently large, the dislocation travels through various Peierls wells (for 
example OA and BC are in different Peierls wells, as are I J and KL in the unloading branch). 
It is clear that the r app vs. d curve obtained will depend on the simulation cell size. On the 
other hand, if the boundary term is subtracted off using eq. (0), the results should be found 
to be cell size independent. Also the lattice resistance function computed from various Peierls 
wells must agree (since the lattice resistance is periodic), in addition to being simulation cell size 
independent. 

To examine the validity of our ideas concerning the boundary force, we have examined the 
response of screw dislocations to an applied stress in aluminum and titanium. The point of these 
calculations is to revisit the idea of lattice resistance of dislocations from an atomistic perspective 
while approximately accounting for the finite simulation cell size. 

Our first test case was the simulation of a | (110) type screw dislocation in fee aluminum 
which was performed using the embedded atom potentials developed by Ercolessi and Adams 
(1993). For details on the embedded atom method, see Daw and Baskes (1984). Fig. 4 shows 
the relationship between r app and d (compare with fig. 3). It is clear that this curve is simulation 
cell-size dependent. The lattice resistance function computed using eq. fl3~?|) is given in fig. 5. 
In these calculations the value of the boundary force coefficient was taken to be 6.3. It is seen 
that the computed lattice resistance function is simulation cell-size independent, with a resulting 
Peierls stress of 0.00068/i. 

Our additional test was made on hep titanium. In this case, embedded atom potentials 
developed by Igarashi, Khantha and Vitek (1991) were used. The screw dislocations simulated 
had a Burgers vector of y(1120) on the basal plane. The lattice resistance function is plotted 
in fig. 6. The anisotropy ratio k for Titanium is 1.06 and the boundary force coefficient A for 
this value of k is 6.2. This value for A was found to be too high; and when A was set to 5.5 the 
lattice resistance obtained from different Peierls wells turned out to be consistent. The value of 
the Peierls stress is about 0.00034/i. 

It must be noted that in plotting the lattice resistance functions in figs. 5 and 6 using eq. (ft7|), 
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the negative sign appearing in the equation has been dropped. Thus the resistance function must 
be interpreted as the external stress necessary to overcome the lattice resistance. In both the 
cases discussed above, the lattice resistance functions from different Peierls wells are plotted in 
the same graph. It is thus evident that the computed lattice resistance functions are periodic. It 
is interesting to note that in all the cases, the lattice resistance functions are multi-welled. The 
screw dislocations in all the materials discussed above undergo a splitting reaction into Shockley 
partials and thus the lattice resistance curves are for these extended dislocations. The effect 
of this splitting on the boundary force is discussed in Appendix I. The value of A based on 
isotropic elasticity works well for the case of Al while in the case of Ti the value of A had to be 
set to 5.5 (as opposed to 6.2 given by the analytic calculation) to obtain consistent values from 
different Peierls wells (i.e., for the resistance function to be periodic). 

4 Bow out of a pinned dislocation segment 

The simulation of the bow out of a pinned dislocation segment is an important step towards 
understanding the mechanics of Frank-Read sources. Work in this direction has hitherto been 
restricted to continuum formulations. Foreman (1967) applied Brown's self-stress (Brown 1964) 
approach to simulate dislocation bow out. Assuming that the dislocation behaves like a flexible 
string with a line tension T, a simpler theory is obtained (Nabarro 1967, de Wit and Koehler 
1959). The latter theory has the advantage of simplicity and we adopt it for the interpretation 
of the results of the atomistic simulations. The goals of these simulation are i) to examine the 
validity of the concept of line tension from an atomistic viewpoint and to compute its value 
from "first principles" and ii) to compare the shape of the dislocation bow out predicted by the 
continuum theory with that obtained from atomistics. 

The atomistic simulation of bow out of a pinned segment under an applied stress is achieved 
as follows. A straight screw dislocation is placed along the axis of a cylindrical cell of atoms of 
length L (the axis of the cylinder is the x 2 -axis in fig. 1), with the radius of the dynamic region 
M set to R. A strain corresponding to the applied stress T app is imposed on this configuration 
and the energy is minimised with respect to the atomic coordinates in the dynamic region. In 
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the energy minimisation step, in addition to the atoms in region F, the atoms on the x 2 = and 
x 2 = L planes are held fixed imitating the pinning effect on the dislocation. 

To obtain the equilibrium shape of the bow out as predicted by continuum theory, we consider 
a screw dislocation on the x 2 axis pinned at x 2 = and x 2 = L, i.e., the length of the pinned 
segment is L. Let the equilibrium shape of the dislocation be given by the function x\ = f(x 2 ) 
(cf. fig. 7). The equilibrium shape is given by the function that minimises the functional 



/(/) = £ (r/l + /* - r app bf + ^ dx 2 , (38) 

/' is the derivative of / with respect to x 2 and A = Afib 2 / (4ir 2 R 2 ), where A is the boundary force 
coefficient. The quadratic term accounts for the additional energy due to the fixed boundary 
and is absent from a standard continuum formulation. It is assumed in this formulation that the 
applied stress is larger that the Peierls stress so that the lattice resistance (misfit energy) can 
be neglected. In addition, the applied stress r app is also assumed to be small compared to the 
shear modulus /x of the material. The function / satisfies the Euler-Lagrange equation of the 
variational principle eq. (BSD, i.e., 



d I f \ . ,T app b A 



dx9 \ . I\ _|_ f 2 I / T 



7zf) = 0, /(0) = f(L) = 0, (39) 



the solution of which is (assuming \ f'\ <^ 1 since r app <^ fi) 

f( v r app b ( sinh y/f(L-x 2 ) + sinhy/fx 2 \ 

f{x 2 ) = —— 1 = v - (40) 

A V sinh^L J 

which expresses the fact that under the influence of the boundary stress the shape of the bow 
out is given by a catenary. Noting that the maximum bow out d occurs at x = L/2, we find that 



d = /(L/2) = ^ 



sechl^JA 



(41) 



The following strategy is adopted in computing the line tension for the results of the atomistic 
simulations. The maximum bow out d is "measured" from the results of an atomistic simulation. 



If d is known, the line-tension T can be obtained from eq. ([y]) as 



T=~. lt f AJ (12) 

4 (sech -1 (1 - 
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If the effect of the boundary is not explicitly accounted for (A — > 0), then d is related to T using 

- sech 



d = lim Tappb 



A->0 A 



L X 
2VT 



t hT 2 

app (43) 



8T 



and one may obtain an uncorrected estimate of the line tension as 

t hT 2 

u ~ 8d ' 1 ] 

These results were used to perform atomistic simulations of dislocation bow out for f (110) 
type screw dislocation in Al using the EAM potentials developed by Ercolessi and Adams (1993). 
The values of the line tension computed for various lengths of the dislocation at various levels 
of applied stresses are tabulated in tables 1, 2 and 3 (in all these cases the value of R is 40. OA). 
The values of line tension computed using eq. (fj2j), which incorporate the boundary stress, are 
length independent. Also, the values do not depend strongly on the values of the applied stress. 
Table 4 shows the computed value of the line tension using a larger simulation cell (R = 57. OA). 
It is seen that the computed values agree well with those computed using the smaller cell. In 
addition, it is seen that the values of line tension computed from eq. (|42|) and those computed 
from eq. (|44] ) agree when the cell size becomes larger and when the length of the dislocation 
is made smaller, as is evident from tables 1 and 4. At larger lengths of the dislocation the 
bow out d becomes larger (for a given value of T app ) and the quadratic approximation eq. fl3"3|) 
for the additional energy is an underestimate and hence the boundary force is underestimated. 
Consequently, the line tension is overestimated, as seen in table 3. A second interesting trend 
is seen in the estimates for the line tension in the case where the boundary contribution is not 
accounted for. Here as the dislocation length increases, so too does the apparent line tension. 
This effect may be attributed to the boundary force which mimics the effects of a higher line 
tension preventing the development of a full bow out. 

It is also useful to compare the shape of the bow out as predicted by the continuum formulation 
(eq. (|40D) and that obtained from atomistics. It is seen that the shape predicted by the continuum 
theory which accounts for the boundary force is in close agreement with that obtained from 
atomistics (cf. fig. 8). Thus the continuum formulation (eq. (|38"D) provides a satisfactory model 
for the study of dislocation bow out. In addition the concept of line tension is found to be 
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meaningful from an atomistic perspective allowing for a determination of the line tension from 
"first principles". The value of line tension is found to be about 0.42/i6 2 which is close to the 
conventional approximation of 0.5/Z& 2 (Nabarro 1967). It must be noted that this dislocation 
undergoes a splitting reaction into Shockley partials as described before and thus the computed 
line tension is for the extended dislocation. It is noted that even with the boundary force 
correction, the computed values of the line tension are slightly different for various levels of 
applied stress. This is believed to arise from the neglect of the lattice resistance term. 

Fig. 9 shows a plot of [w 2 ] as a function of the position on the slip plane extracted from the 
results of an atomistic calculation using the technique given in Miller and Phillips (1996). The 
value of d is found by sectioning this surface at a value of one half of the Burgers vector. The 
positions of atoms on the atomic planes that sandwich the slip plane is are plotted in fig. 10. A 
careful study of this figure brings out the discrete nature of the dislocation bow out and the core 
structure of a curved dislocation. 

5 Conclusion 

In this paper, we have described a method for accounting for the additional boundary stresses 
that act on a dislocation due to the use of a finite sized cell in an atomistic simulation. An 
explicit formula for the boundary stress is derived for the case of the screw dislocation. The 
formula was derived based on the Peierls-Nabarro model and a linear elastic reckoning of the 
elastic energy due to an inconsistent boundary. The boundary stress is proportional to the inverse 
square of the size of the simulation cell. This formulation was then used to obtain the lattice 
resistance function for screw dislocations in several close packed metals. In addition, a simple 
3D configuration for the bow out of a screw dislocation was also studied and the line tension of 
the dislocation was computed for the case of aluminum. 

The explicit accounting of the boundary stress is essential to obtain meaningful quantitative 
values from the results of an atomistic simulation. This is evident from the simulation of the 
dislocation bow out that is presented here. The effects of the far-field boundary conditions may 
be overcome with the use of large simulation cells when empirical or semi-empirical potentials are 
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used in the simulation, however it proves to be a serious difficulty in calculations using density 
functional theory. Thus this method might be especially useful in those simulations where there 
are severe restrictions on the size of the simulation cell. 



Appendix : Effect of dislocation splitting on the boundary 
force 

It was noted in section 3 that the screw dislocations in aluminum and titanium undergo a splitting 
reaction. Here we compute the boundary force when the dislocation is split assuming a splitting 
distance of 2s ( s = 4.0A for aluminum and s = 6.1 A for titanium). It is also assumed that 
this distance does not change as the dislocation moves under the influence of the applied stress 
- a premise that appears to be in good agreement with observations. Under these assumptions 
one can compute the boundary stress to be 

T ^ d ) = A7[ 2 R 2 d ( 45 ) 

for d <C R, where c = s/R. This relation can be derived using equations eqs. ||, |9| and [LC] where 
Up is taken to be the sum of the displacement fields due to two screw dislocations with Burgers 
vector b/2 located at x\ = d — s and x\ — d + s. In the derivation the edge part of the Burgers 
vector of the Shockley partials is neglected, and thus eq. (flop is approximate. It is seen that the 
boundary stress coefficient depends on the splitting ratio c in addition to the anisotropy ratio 
k. A plot of the boundary force coefficient as a function of the splitting ratio is given in fig. 11. 
The value of A is seen to be unaffected when the splitting is small (c < 0.1). This condition is 
satisfied in the cases considered here and the effect of splitting is thus small. 
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Figure Captions 



Figure 1: Geometry of the simulation cell illustrating the implementation of a far field boundary 
condition by freezing atoms in positions dictated by linear elasticity theory. 

Figure 2: Dimensionless stress F(a, k) as a function of dimensionless displacement of dislocation 
line, shown for different values of the anisotropy ratios. 

Figure 3: Schematic of r app vs. d curve. (d a is the repeat lattice distance) 

Figure 4: r app vs. d as obtained from atomistic simulations for Al (compare with fig. 3) showing 
size dependent results in the absence of boundary correction. 

Figure 5: Lattice resistance function for screw dislocation in Al accounting for boundary stress 
and showing system size independence. Dashed lines are provided as a schematic indication 
of the parts of the resistance curve that are inaccessible using this method. 

Figure 6: Lattice resistance function for screw dislocation in Ti. Dashed lines are provided as 
a schematic indication of the parts of the resistance curve that are inaccessible using this 
method. 

Figure 7: Geometry of dislocation bow out. 

Figure 8: Comparison of atomistic and continuum bow out in Al. R = 40. OA, L = 211. 2A, 
r app = 0.0025//. 

Figure 9: Plot of [u 2 ] on the slip plane. R = 40.0A, L = 211.2A, r app = 0.01//. 

Figure 10: Positions of atoms on atomic planes that sandwich the slip plane. R = 40. OA, 
L = 211.2A, r app = 0.01//. 

Figure 11: Variation of the boundary stress coefficient A with normalised splitting distance c. 
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d(A) 


T u / fib 2 (without boundary effects) 


T / fib 2 (with boundary effects) 


0.0025 


0.83 


0.43 


0.40 


0.0050 


1.59 


0.45 


0.42 


0.0100 


3.40 


0.42 


0.39 



Table 1: Bow out and line tension for screw dislocation in Al. (L = 57.2A, R =40.0A) 





d (A) 


T u j fib 2 (without boundary effects) 


T / fib 2 (with boundary effects) 


0.0025 


2.43 


0.53 


0.41 


0.0050 


4.66 


0.55 


0.43 


0.0100 


9.11 


0.56 


0.44 



Table 2: Bow out and line tension for screw dislocation in Al. (L = 108. 5 A, R =40. OA) 



T~app/ fJ> 


d (A) 


T u j fib 2 (without boundary effects) 


T j fib 2 (with boundary effects) 


0.0025 


5.19 


0.94 


0.46 


0.0050 


10.24 


0.95 


0.48 


0.0100 


18.74 


1.04 


0.57 



Table 3: Bow out and Line Tension for Screw Dislocation in Al. (L = 211.2A, R =40.0A) 



7~appl 


d (A) 


T u j fib 2 (without boundary effects) 


T j fib 2 (with boundary effects) 


0.0025 


0.88 


0.42 


0.40 


0.0050 


1.62 


0.46 


0.44 


0.0100 


3.47 


0.43 


0.42 



Table 4: Bow out and line tension for screw dislocation in Al. (L = 58.1 A, R =57.0A) 
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Shenoy and Phillips, Figure 1 




Shenoy and Phillips, Figure 2 




Shenoy and Phillips, Figure 3 




Shenoy and Phillips, Figure 5 




Shenoy and Phillips, Figure 6 




Shenoy and Phillips, Figure 7 



Shenoy and Phillips, Figure 8 




Shenoy and Phillips, Figure 9 
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Shenoy and Phillips, Figure 10 




Shenoy and Phillips, Figure 11 



